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Naval  Ocean  Research  and  Development  Activity  Code  323 
National  Space  Technology  Laboratories.  Mississippi 

Abstract  A  skillful  ocean  fo^a^mg  capability  will  de¬ 
pend  on  W  adequate  data  inputs  fSV-adequate^computing 
power,  and  'f^f'properly  designed  and  adequately  validated 
ocean  models  for  data  assimilation  and  forecasting.  Once  these 
conditions  are  satisfied, -^hich  is  feasible  with  existing  tech¬ 
nology,  forecasts  of  meandering  currents  and  eddies  up  to  a 
few  months  appear  to  bea  reasonable  expectation.  Simulation 
studies  suggest  that  tmprediction  of  meandering  currents  anct^  ^  - 
eddies  requires  high  horizontal  resolution UffOkm), but  only> 

-  low  vertical  resolution.  High  vertical  resolution  can  be  ob*  v 
tained  by  using  the  circulation  model  to  provide  horizontal  and 
vertical  advection  to  a  grid  of  one-dimensional  mixed-layer 
models.  Circulation  models  on  subdomains  of  major  ocean  ba- 
sins  with  extensive  open  ocean  boundaries  require  accurate  /r 
boundary  information  at  all  levels  in  the  verticid  and  for  the  / 
duration  of  forecasts  longer  than  a  few  days.  Basin-scale 
models  appear  to  be  the  most  promising  source  for  this  infor-  / 
mhtion.  To  provide  useful  boundary  conditions,  the  basin-scale 
models  must  also  resolve  the  meandering  currents  and  eddies.  / 

Class  7  computers  (-^Tgigaflopand  32—128  million  words), 
expected  on  the  market  in  theUfiddlc  and  hKM9t0s,  are  re¬ 
quired  for  eddy-resolving  models  of  nqjor  ocean  basins. Cur¬ 
rently,  satellites  and  atmospheric  models  provide  the  only 
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prospects  for  oceanographic  data  and  forcing  functions  with 
global  coverage  and  resolution  adequate  for  ocean  circulation 
prediction  models.  For  this  purpose  the  most  useful  fields  po¬ 
tentially  available  from  satellites  are  current-related  variations 
in  the  sea-surface  elevation  derived  from  altimeter  data  and 
scatterometer-derived  wind  stress.  Choosing  an  appropriate  al¬ 
timeter  track  pattern  is  critical.  Successful  ocean  prediction 
appears  feasible  without  extensive  subsurface  data  acquisition 
\by  using  the  circulation  model  to  convert  the  potentially  well- 
observed  fields  at  the  surface  into  subsurface  information.  This 
looks  more  promising  than  attempting  to  assimilate  the  ex¬ 
tremely  sparse  subsurface  data.  Where  adequate  subsurface 
data  are  available  (nowhere  at  present),  they  should,  of  course, 
be  used.  Subsurface  data  and  major  field  programs  are  vital 
for  local  and  regional  subsurface  validation  of  the  forecast 
models,  but  altimeter  data  are  essential  for  global  prediction 
and  representation  of  individual  current  meanders  and  eddies. 
Otherwise,  in  most  regions  it  would  be  necessary  to  rely  on 
the  simulation  skill  of  the  model  and  the  forcing  functions, 
which  are  only  weakly  coupled  to  many  of  these  oceanographic 
features.  Within  the  next  decade  the  appropriate  satellite  data. 
Class  7  computers,  and  eddy-resolving  basin  or  global-scale 
models  should  become  available,  and  should  be  used  to  form 
the  heart  of  a  global  ocean  prediction  system.  Part  of  this  paper 
outlines  a  strategy  for  global  ocean  prediction  based  on  these 
elements  and  on  plans  at  the  Naval  Ocean  Research  and  De¬ 
velopment  Activity. 


latrododfcNi 

The  next  decade  should  be  to  numerical  prediction  of  ocean  cir¬ 
culation  what  the  1950s  and  1960s  were  to  numerical  weather 
prediction:  a  time  when  the  essential  dements  came  together  to 
permit  major  advances.  These  elements  me  (1)  adequate  data 
input,  (2)  adequate  computing  power,  and  (3)  properly  designed 
and  adequately  validated  models  for  both  data  assimilation  and 
forecasting.  Ten  years  ago  none  of  these  requirements  were  sat¬ 
isfied  for  many  aspects  of  ocean  forecasting.  Ilsandsring  cur¬ 
rents  and  eddies  provide  a  notable  example:  dear  prospects  for 
adequate  data  with  high  density  and  global  coverage  from  sat- 
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ellite-bome  instruments  have  become  evident  only  in  the  last 
few  years;  Class  7  computers  with  a  sustainable  speed  of  - 1 
gigaflop  and  32  to  128  millions  words  of  memory  are  required 
for  eddy-resolving  forecasts  in  major  ocean  basins;  and  notable 
differences  from  meteorology  are  required  in  model  design  and 
in  data  assimilation. 

The  U.S.  Navy  has  determined  a  requirement  for  an  ocean 
prediction  capability,  and  has  established  a  research  program  to 
meet  the  requirement.  This  paper  addresses  the  requirements  for 
data  and  computing  power  and  the  issues  of  model  design  and 
data  assimilation,  primarily  in  the  context  of  the  Navy  effort  at 
the  Naval  Ocean  Research  and  Development  Activity 
(NORDA);  it  is,  however,  an  overview  and  not  a  report  on  any 
specific  research  project.  It  should  be  stressed  that  there  are 
many  types  of  ocean  prediction,  aimed  at  different  aspects  of 
the  motion,  internal  structure,  and  surface  of  the  ocean.  This 
paper  is  focused  on  one  important  aspect  where  altimeter  data 
can  play  a  central  role.  In  particular,  the  role  that  satellite  alti¬ 
meter  data  can  play  will  be  demonstrated  for  detailed  surface 
and  subsurface  forecasts  of  meandering  currents,  eddies,  and 
frontal  locations,  forecasts  that  may  extend  up  to  several 
months. 

Different  Types  of  Ocean  Prediction 

Before  concentrating  on  a  particular  type  of  ocean  prediction, 
we  will  briefly  survey  the  range  of  ocean  phenomena  for  which 
forecasts  are  desirable.  If  we  focus  on  a  particular  region  of  the 
ocean,  then  the  evolution  of  the  motion,  internal  structure,  and 
surface  from  time  f0  to  time  it  depends  oa  (1)  the  initial  state  at 
to,  (2)  the  atmospheric,  astronomical,  and  geological  forcing 
functions,  (3)  flow  through  any  open  boundary  segments  of  the 
region,  (4)  the  bottom  topography  and  coastline  geometry,  and 
(5)  the  physical  and  dynamic  properties  of  the  ocean  fluid.  The 
ocean  exhibits  a  wide  variety  of  phenomena  on  a  broad  spectrum 
of  time  and  space  scales  that  differ  greatly  in  their  manner  and 
degree  of  dependence  on  each  of  these  factors.  Thus  in  predicting 
different  phenomena  we  anticipate  (1)  different  degrees  of  ac- 
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curacy  in  the  initial  state,  (2)  different  time  scales  for  accurate 
forecasts,  (3)  different  modeling  strategies,  (4)  different  data  re¬ 
quirements,  and  (5)  different  data  acquisition  and  sampling  strat¬ 
egies.  Also,  note  that  these  five  aspects  can  depend  on  inter¬ 
actions  between  phenomena,  time  and  space  scales,  and,  of 
course,  the  adequacy  of  the  models  and  the  data  input.  Fortu¬ 
nately,  many  phenomena  and  scales  are  sufficiently  decoupled 
that  it  is  possible  to  design  forecast  models  that  can  simulate 
and  forecast  certain  phenomena  while  suppressing  or  parame¬ 
terizing  the  statistical  effects  of  the  remainder. 

Table  1  shows  some  useful  categories  of  oceanic  response  to 
atmospheric  forcing  that  cover  many  phenomena  of  potential 
interest  in  ocean  forecasting.  It  is  not  all-encompassing  (e.g., 
tides,  tsunamis,  individual  salt  fingers,  and  convection  cells  are 
omitted),  but  is  intended  to  help  place  the  multifaceted  problem 
of  ocean  forecasting  into  perspective,  identify  phenomena  where 
predictive  skill  is  feasible,  clarify  the  potential  roles  of  altimeter 
data,  and  define  the  focus  of  the  paper,  which  is  primarily  on 
Class  2  as  defined  in  the  table. 

“Nowcasting”  is  important  in  all  three  classes  listed  in  Table 
1.  To  nowcast,  the  models  are  integrated  forward  in  time  while 
being  driven  by  new  atmospheric  forcing  functions  and  assimi¬ 
lating  appropriate  new  oceanic  data  as  they  become  available. 
This  allows  the  models  to  fill  in  temporal  gaps  in  the  data  by 
using  their  predictive  skill,  convert  better-observed  surface 
fields  into  subsurface  structure,  and  convert  better-observed  at¬ 
mospheric  forcing  functions  into  oceanographic  information. 

Up  until  now,  ocean  forecasting  efforts  have  concentrated  pri¬ 
marily  on  tides  and  certain  phenomena  in  Class  1,  a  class  which 
to  some  extent  includes  the  cumulative  effects  of  fine  scale  phe¬ 
nomena  predicted  by  empirical  or  semtemptrical  parameteriza- 
tkrns.  The  U.S.  Navy  has  recently  initiated  hemispheric  nu¬ 
merical  forecasts  of  the  upper  mixed  layer  at  the  Fleet  Numerical 
Oceanography  Center  (Clancy  and  Martin,  1981;  Clancy  and  Pol¬ 
iak,  1983).  It  has  also  installed  a  hemispheric  surface  wave  fore¬ 
cast  model  (Pierson,  1982).  Among  ocean  phenomena,  numerical 
storm  surge  prediction  has  a  relatively  long  history  (Welander, 
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1961 ;  Jelesnianski,  1967;  Crawford,  1979).  Coastal  upwelling  has 
been  the  subject  of  a  few  hindcasting  studies  (Hebum,  1980). 

Although  a  satellite  altimeter  can  measure  surface  wind  speed 
and  significant  wave  height  in  addition  to  sea-surface  elevation 
(Fedor  and  Brown,  1962),  the  value  of  a  single  altimeter  with 
only  a  nadir  beam  is  severely  degraded  when  applied  to  predic¬ 
tion  of  phenomena  in  Class  1.  This  is  due  to  the  inverse  relation 
between  spatial  and  temporal  resolution,  which  prevents  ade¬ 
quate  resolution  for  these  phenomena  in  both  space  and  time 
simultaneously.  This  is  also  true  for  tsunamis,  which  have  time 
scales  of  a  few  hours  (Neumann  and  Pierson,  1966).  However, 
the  altimeter  can  be  useful  in  constructing  climatologies  of  wind 
speed  and  significant  wave  height  (Chehon  et  al.,  1981;  Wentz 
et  al.,  1962;  Mognard  et  al.,  1963),  in  validation  and  tuning  of 
certain  Class  1  forecast  models,  and  in  some  research  problems 
involving  these  quantities.  Because  tides  are  repetitive  with 
known  periods,  the  altimeter  should  be  important  in  refining  the 
global  knowledge  of  tidal  phase  and  amplitude  (Cartwright  and 
Alcock,  1961;  Brown  and  Hutchinson,  1961;  Diamante  and  Nee, 
1961;  Mazzega,  1963). 

The  output  from  atmospheric  prediction  models  and  satellite- 
borne  instruments  measuring  surface  wind  speed  and  direction 
(scatterometers)  or  even  wind  speed  done  (scanning  multifre¬ 
quency  microwave  radiometers  [SMMRs])  shows  great  promise 
in  facilitating  predictions  in  Class  1.  Scatterometers  and  micro¬ 
wave  radiometers  can  provide  adequate  coverage,  accuracy .  and 
spatial  and  temporal  resolution,  except  for  coverage  near  coastal 
boundaries  (Jones  et  al.,  1962;  Wentz  et  al.,  1962;  Lipes,  1982; 
Satellite  Surface  Stress  Working  Group,  1962;  Mueller,  1962). 
For  mixed-layer  forecasting,  amasurements  of  sea-surface  tem¬ 
perature  (SST)  (multichannel  infrared  and  microwave  radiome¬ 
ters)  would  also  be  of  value.  However,  SST  measurements  re¬ 
main  plagued  by  cloud  contamination  (iafared  radiometers)  or 
relatively  poor  accuracy  and  resolution  (SMMR,  123  km  at  6.6 
GHz  with  1-m  diameter  antenna;  25  km  at  6.6  GHz  with  the  4- 
m  diameter  parabolic  dish  antsnns  proposed  for  die  canceled 
National  Oceanic  Satellite  System,  NOSS)  (Mueller,  1962).  Res- 
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olution  of  25  km  is  adequate  for  most  mesoscaie  phenomena. 
For  both  wind  speed  and  SST  the  SMMR  on  SEASAT  (1-m 
diameter  antenna)  yielded  poor  results  in  areas  of  sun  glint  and 
within  600  km  of  land  (Lipes,  1982).  Some  work  has  also  been 
done  on  satellite  measurement  of  latent  heat  fluxes  (Liu,  1984) 
and  incident  solar  radiation  at  the  surface  (Gautier,  1981),  the 
two  heat  fluxes  that  are  usually  the  largest  at  the  air-sea  interface. 

The  prediction  of  mesoscaie  (50  to  500  km)  eddies,  meandering 
currents,  and  frontal  positions  (phenomena  in  Class  2)  is  the 
subject  of  the  remaining  sections  of  this  paper.  The  prediction 
of  this  class  of  phenomena  is  one  area  where  altimeter  data  can 
play  a  central  role.  This  is  because  (1)  individual  mesoscaie  ed¬ 
dies  and  current  meanders  are  often  not  driven  directly  by  the 
wind  or  by  any  other  external  forcing  function,  and  (2)  eddies 
can  persist  for  more  than  a  year  after  their  initial  generation  (Lai 
and  Richardson,  1977).  Thus,  oceanic  data  are  crucial  for  reliable 
prediction  of  individual  eddies  and  current  meanders.  With  suf¬ 
ficient  oceanic  data,  we  anticipate  that  the  prediction  of  these 
features  can  be  treated  as  an  initial  value  problem  in  which  the 
future  forcing  functions  are  representative  but  not  accurate  on 
time  scales  greater  than  a  few  days.  Without  oceanic  data  input, 
simulation  can  be  used  to  predict  the  statistical  properties  of 
features  and  ensembles,  but  usually  not  the  evolution  and  move¬ 
ment  of  the  individual  features. 

Up  until  now  there  has  been  little  work  on  the  prediction  of 
Class  2  phenomena,  except  for  the  work  on  the  prediction  of 
mesoscaie  eddies,  using  limited-area  models,  by  A.  R.  Robin¬ 
son's  group  at  Harvard  (Robinson  and  H aid vogel,  1980;  Robin¬ 
son  mtd  Tu,  1982;  Miller  and  Robinson,  1984).  However,  there 
is  a  substantial  body  of  literature  on  the  simulation  and  ocean 
dynamics  of  Class  2  phenomena  that  is  very  useful  in  the  design 
of  ocean  circulation  prediction  models  (Holland  and  Lin,  1975; 
Rhines,  1977;  Robinson  et  al.,  1977;  Semtner  and  Mintz,  1977; 
Holland,  1978, 1982;  Semtner  and  Holland,  1978;  McWilliams  et 
al.,  1978;  McWilliams  and  Flier!,  1979;  Cox,  1979;  Holland  and 
Rhines,  1980;  Hurlburt  and  Thompson,  1980, 1982;  Lin  and  Hurl¬ 
burt,  1981;  Schmitz  and  Holland,  1982;  Heburn  et  al.,  1962). 
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ployed  by  Sarmiento  and  Bryan  (1982)  and  Clancy  and  Poliak 
(1983). 

Thus  far  only  nowcasting  of  Class  3  phenomena  has  been  dis¬ 
cussed.  Forecasting  of  these  phenomena  is  based  on  the  hy¬ 
pothesis  that  although  they  are  a  direct  response  to  the  atmos¬ 
pheric  forcing  functions,  the  response  is  slow  enough  to  be 
insensitive  to  forecast  errors  in  the  daily  fluctuations  of  the  at¬ 
mosphere.  This  insensitivity  should  permit  forecasts  on  time 
scales  greater  than  the  few  days  possible  for  atmospheric  pre¬ 
diction. 

Haney  (1980)  has  performed  non-eddy-resolving  nowcasting 
and  forecasting  studies  of  large-scale  ocean  anomalies,  using  a 
closed  rectangular  domain  covering  a  major  portion  of  the  North 
Pacific,  a  region  adequately  sampled  for  large-scale  anomalies 
at  the  time  by  ships  of  opportunity  (White  mid  Bernstein,  1979). 
The  model  included  simple  mixed-layer  physics  and  ten  levels 
in  the  vertical  with  six  in  the  upper  262  m.  The  model  was  spun 
up  for  230  years  by  using  climatological  forcing.  Then  four  ex¬ 
periments  were  performed.  In  the  first  an  observed  thermal 
anomaly  was  added  to  the  climatoiogically  driven  state  and  the 
model  was  integrated  for  four  additional  months  by  using  cli¬ 
matological  forcing.  The  four  month  forecast  was  poor,  but 
slightly  better  than  a  forecast  of  persistence.  Next,  three  nowcast 
experiments  (in  the  sen*  described  earlier)  were  performed. 
Each  time  the  model  teas  integrated  over  four  months  by  using 
observed  forcing  fonctioni,  which  had  been  averaged  for  a 
month.  The  options  wore  observation!  or  dmmtotogy  for  the 
winds,  the  hem  Bunas,  and  the  initial  stale,  la  the  first  nowcast 
experiment,  observations  were  used  only  for  the  winds,  in  the 
second  only  for  the  hem  Duxes,  end  in  dm  third  for  afi  three 
foncrtoei.  h  el  three  cnees  the  results  were  quite  good  in  the 
upper  NO  m  and  much  better  time  pertinence  or  the  forecast. 
Between  MO  and  300  m  fee  results  seers  not  as  good,  but  better 
than  persistence  and  better  than  the  forecast,  except  m  the  sec¬ 
ond  experiment,  where  the  results  wove  poor,  dearly,  additional 
research  is  required  to  determine  the  neifolaiis  of  the  nowcast 
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Many  of  these  have  simulated  some  features  of  the  ocean  cir¬ 
culation  with  substantial  success  by  using  simple  models,  rec¬ 
tangular  domains,  and  simple  steady  forcing  functions.  These 
successes  are  an  indication  of  potential  predictive  skill,  but  one 
that  is  largely  unverified  at  this  time.  It  is  noteworthy  that  no 
ocean  modeling  study  yet  published  has  included  all  of  the  fol¬ 
lowing  basic  features:  (1)  resolution  sufficient  for  mesoscale  ed¬ 
dies,  (2)  stratification,  (3)  realistic  coastline,  (4)  realistic  bottom 
topography,  (5)  thermodynamics  and  mixed-layer  physics,  (6) 
integration  to  statistical  equilibrium,  and  (7)  a  deep  ocean  basin 
with  dimensions  (L)  large  enough  (L  ~  1000  km  or  more)  for 
planetary  vorticity  advection  to  be  important  (that  is,  L  >  ( VI 
0)'72,  where  V  is  a  characteristic  velocity,  0  is  the  variation  of 
the  Coriolis  parameter /with  latitude  4>,/  «  2u>  sin  4>,  and  u>  is 
the  angular  velocity  of  the  earth's  rotation. 

Phenomena  in  Class  3  are  a  direct  response  to  atmospheric 
forcing,  but  on  much  longer  time  scales  than  those  in  Class  1 . 
Because  they  are  a  much  more  integrated  response  to  the  wind, 
the  requirements  for  temporal  resolution  and  accurate  depiction 
of  daily  atmospheric  fluctuations  are  not  as  stringent  as  in  Class 
1 .  In  Class  3  “nowcasting"  and  forecasting  may  be  feasible  from 
new  atmospheric  forcing  functions  without  new  oceanic  data. 

Class  3  predictions  might  be  accomplished  by  spinning  up  the 
forecast  model  to  the  present,  using  the  best  historical  data  avail¬ 
able  for  the  forcing  functions.  The  spin-up  should  be  long  enough 
to  allow  the  model  to  reach  statistical  equilibrium.  Busaiacchi 
et  a).  (1983)  have  applied  this  technique  to  a  simple  linear  model 
of  the  tropica]  Pacific  Ocean,  with  some  success  on  seasonal  and 
interannual  time  scales.  Preller  and  Hebum  (1983)  have  shown 
that  a  simple  model  of  the  western  Mediterranean  Sea,  driven 
by  mean  monthly  winds  from  May  (1982),  can  reproduce  many 
of  the  mqjor  persistent  features  observed  with  scales  greater  than 
100  km.  The  addition  of  altimeter-derived  sea-surface  elevations 
accurate  on  the  scales  of  mqjor  ocean  basins  as  well  as  on  me- 
soscales  would  reduce  the  burden  on  the  simulation  skill  of  the 
model.  Relaxation  of  the  mode)  simulation  toward  observed  cli¬ 
matology  would  also  reduce  the  burden.  This  is  a  technique  em- 
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Time  and  Space  Scales  for  the  Prediction  of 
Mqjor  Current  Systems  and  Eddies 

The  U.S.  Navy  is  particularly  interested  in  phenomena  that  can 
substantially  affect  acoustic  propagation  on  tactical  space  scales 
up  to  several  hundred  kilometers.  As  a  result  this  paper  will  focus 
on  the  potential  for  forecasting  mqjor  current  systems,  mesoscale 
eddies,  frontal  locations,  and  the  mixed  layer,  all  features  that 
have  been  demonstrated  to  have  an  important  impact  on  acoustic 
propagation  (Baer,  1980;  Beckerle  et  al.,  1980;  Nysen  et  at.,  1978; 
Lawrence,  1983;  Itzikowitz  et  al.,  1983). 

A  schematic  diagram  from  the  report  of  the  TOPEX  Science 
Working  Group  (1981)  is  quite  instructive.  Figure  1  plots  the 


FIGURE  1.  Schematic  diagram  of  the  time  scale  versus  the  space  scale  for 
oceanic  variability  away  from  boundary  layers  and  contoured  in  arbitrary  units 
(from  the  TOPEX  Science  Working  Group,  1961). 
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amplitude  of  oceanic  variability  versus  time  scale  and  horizontal 
length  scale.  It  is  representative  of  the  open  ocean  tens  of  meters 
from  the  surface  or  bottom  and  a  few  hundred  kilometers  from 
the  equator  or  coastal  boundaries.  The  peak  in  variability  around 
sixty  days  and  100  km  associated  with  Class  2  phenomena  implies 
the  need  for  about  10- km  grid  resolution  in  the  corresponding 
numerical  ocean  circulation  forecast  models  (Grammeltvedt, 
1969,  Haidvogel  et  al.,  1980).  A  similar  peak  for  the  atmosphere 
occurs  at  about  six  days  and  3000  km  (Vinnichenko,  1970;  Kao 
and  Lee,  1977),  indicating  that  this  resolution  is  dynamically 
comparable  to  300  km  in  atmospheric  models.  Atmospheric  pre¬ 
dictability  has  received  substantial  attention  (e.g.,  Lorenz,  1982). 
Usable  predictive  skill  of  about  one  week  is  indicated  for  at¬ 
mospheric  systems,  systems  dominated  by  major  cyclones,  an¬ 
ticyclones,  and  waves  on  the  jet  stream.  Comparable  predicta¬ 
bility  studies  have  not  been  made  for  the  ocean.  However,  as 
noted  by  Haidvogel  and  Holloway  (1984),  predictability  studies 
using  simple  quasigeostrophic  models  can  be  applied  to  either 
the  ocean  ( o )  or  the  atmosphere  (a)  by  appropriate  scaling,  and 
a  reasonable  choice  of  velocity  scales  (V0  ~  10  cm/s,  Va  ~  10 
m/s)  and  length  scales  (La  ~  100  km,  La  ~  1000  km)  yields  an 
oceanic  time  scale  ( L0IV0 ),  O(10x)  the  atmospheric  one.  This 
is  consistent  with  McWilliams  and  Shen  (1980),  who  suggest 
2itL0IV0  =*  60  days  as  a  crude  estimate  of  the  predictability  time 
scale  for  a  field  of  mesoscale  ocean  eddies.  The  quasigeostrophic 
scaling  and  the  observed  six-  versus  sixty-day  peak  for  atmos¬ 
pheric  versus  open  ocean  temporal  variability  at  least  hint  that 
oceanic  forecast  skill  for  Class  2  phenomena  may  extend  about 
ten  times  longer  than  atmospheric  predictive  skill. 

Quite  different  results  have  been  obtained  by  Miller  and  Rob¬ 
inson  (1984).  They  have  used  a  six-layer  quasigeostrophic 
model  with  a  500-km  by  500-km  open  domain  to  study  the  pre¬ 
dictability  of  eddies  in  the  region  of  the  POLYMODE  field  pro¬ 
gram  in  the  Atlantic  southwest  of  Bermuda.  This  is  a  relatively 
benign  region  with  no  mqjor  current  systems  passing  through  it. 
Instead  it  is  populated  by  mesoscale  eddies  that  are  consistent 
with  Figure  1  and  that  have  flow  speeds  around  20  cm/s  and 
westward  propagation  speeds  of  about  6  cm/s.  The  interior  and 
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the  open  boundaries  of  the  forecast  model  were  initialized  at  all 
levels  in  the  vertical  by  error-free  simulated  data  taken  from  an 
integration  of  the  same  model  on  a  larger  domain,  an  integration 
also  used  for  forecast  verification.  Atmospheric  forcing  was 
omitted.  When  the  open  boundary  conditions  were  held  fixed  in 
time.  Miller  and  Robinson  found  that  the  forecast  was  superior 
to  persistence  for  the  interior  field  for  only  about  two  days  (R. 
N.  Miller,  Tulane  University,  New  Orleans,  personal  commu¬ 
nication,  1983)  and  that  errors  reached  the  30%  level  in  about 
four  days.  However,  errors  remained  low  for  the  three-month 
duration  of  the  experiments,  if  the  open  boundary  conditions 
were  updated  frequently  enough  (at  least  every  four  days  in  their 
case).  Carter  and  Robinson  (1981)  have  performed  a  space-time 
objective  analysis  of  the  depth  of  the  15°C  isotherm  (mean  depth 
=*  600  m)  using  POLYMODE  expendable  bathythermographs 
(XBTs).  The  analysis  was  dominated  by  mesoscale  features  and 
showed  that  a  forecast  of  persistence  would  yield  a  linear  cor¬ 
relation  of  .6  after  ten  days,  .3  after  twenty  days,  and  - .  1  after 
forty  days. 

The  results  in  the  two  preceding  paragraphs  are  useful  in  de¬ 
signing  a  system  for  oceanic  prediction.  The  results  from  Miller 
and  Robinson  (1984)  do  not  apply  directly  to  enclosed  or  semi- 
enclosed  domains  or  to  much  larger  domains  with  open  bound¬ 
aries,  because  in  their  case  the  open  boundary  conditions  are 
much  more  important  than  the  initial  conditions  in  determining 
the  time  scales  for  predictability.  However,  their  results  do  em¬ 
phasize  that  great  care  must  be  exercised  with  open  boundary 
segments  or  adjacent  data-poor  regions  (McWilliams  and  Shen, 
1980)  because  of  their  potential  for  extreme  impact.  Research 
on  the  methodology  and  effects  of  open  boundaries  must  play 
an  important  role  in  ocean  prediction  studies,  because  open  or 
artificially  closed  boundary  segments  are  a  necessary  feature  of 
any  ocean  circulation  prediction  model  that  does  not  cover  the 
entire  global  ocean. 

It  should  be  noted  that  the  Naval  tactical  scale  and  the  hori¬ 
zontal  scale  of  maximum  temporal  variability  strongly  overlap. 
Furthermore,  the  energy  of  the  time-varying  circulation  exceeds 
the  energy  of  the  mean  circulation  over  most  of  the  world  ocean 
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(Wyrtki  et al.,  1976).  At  oceanic  scales  greater  than  500  km  the 
temporal  variability  is  relatively  weak  except  for  an  annual  var¬ 
iation  (Figure  1).  Although  large-scale  anomalies  are  significant 
in  some  regions  and  could  be  detected  (White  and  Bernstein, 
1979)  and  at  least  to  some  extent  predicted  (Haney,  1980),  this 
suggests  that  the  coarse  horizontal  grids  used  for  atmospheric 
models  are  of  limited  value  for  ocean  forecasting  except  when 
the  resolved  scales  of  the  atmospheric  forcing  largely  determine 
the  scales  of  the  oceanic  response,  for  example  in  prediction  of 
surface  waves  (Pierson,  1982)  and  changes  in  the  mixed  layer 
(Clancy  and  Poliak,  1983). 

Some  Obstacles  to  the  Development  of  a  Numerical  Ocean 
Prediction  Capability 

A  list  of  desirable  characteristics  for  a  comprehensive  ocean  pre¬ 
diction  model  capable  of  handling  Class  2,  Class  3,  and  certain 
Class  1  phenomena  as  described  in  Table  1  (that  is,  mixed  layers 
and  equatorial  and  some  coastal  phenomena  with  short  time 
scales)  might  include  the  following: 

1.  High  horizontal  resolution,  -10  km. 

2.  Mixed-layer  physics. 

3.  High  vertical  resolution  for  the  mixed-layer  physics  and  for 
the  acoustically  important  temperature  and  salinity  struc¬ 
ture. 

4.  A  limited  domain  (-1000  km  by  1000  km  or  less)  to  cover 
the  tactical  space  scale  of  naval  interest  or  to  provide  other 
-tailored  services  (e.g.,  for  fisheries  or  ice,  oil  spill,  or  other 
pollutant  trajectories)  while  at  the  same  time  minimizing 
the  computation. 

5.  Initialization  from  observational  data  (as  in  meteorology). 

6.  Predicted  forcing  functions  on  the  time  scale  of  their  va¬ 
lidity,  and  representative  but  not  accurate  forcing  functions 
on  longer  time  scales. 

To  a  large  extent  this  list  follows  the  lead  of  meteorology,  a 
science  with  thirty  years  of  experience  in  numerical  prediction. 
However,  a  numerical  ocean  prediction  system  with  all  of  these 
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characteristics  would  encounter  serious  obstacles: 

1.  In  limited-area  models  (500  km  by  500  km).  Miller  and  Rob¬ 
inson  ( 1904)  have  found  that  extensive  open  boundaries  that 
are  poorly  known  pose  the  most  serious  limitation  to  fore¬ 
cast  skill. 

2.  Computing  power  places  a  constraint  on  the  combinations 
of  horizontal  and  vertical  resolution,  domain  size,  and 
model  physics  that  are  feasible  (see  section  on  computing 
requirements). 

3.  At  the  very  least,  forecasts  of  individual  current  meanders 
and  eddies  will  require  adequate  oceanographic  data  as 
input.  This  is  not  provided  on  a  regular  basis  by  the  sources 
for  subsurface  data  currently  available  (see  section  on  data 
requirements). 

Up  until  now,  limitations  in  computing  power  have  forced  se¬ 
rious  compromises  on  eddy-resolving  ocean  models.  For  a  model 
with  10-lun  grid  resolution,  there  has  been  little  choice  but  to 
confront  problems  with  artificial  boundaries,  except  possibly  in 
semienclosed  seas  such  as  the  Mediterranean.  Furthermore, 
forecasting  with  real  data  has  been  limited  to  situations  where 
major  field  programs  such  as  FOLYMODE  have  provided  ad¬ 
equate  oceanographic  data. 

Using  a  barotropic  model  in  a  region  that  is  relatively  simple 
and  homogeneous  and  has  well-known  statistics,  Robinson  and 
Tu  (1982)  have  found  that  an  optimal  estimation  combination  of 
persistent  and  statistically  forecast  boundary  conditions  can  sub¬ 
stantially  reduce  the  error  growth  rate  compared  to  either  alone. 
Embedding  limited-area  models  in  non-eddy-resolving  basin- 
scale  models  is  another  possibility,  but  one  not  very  promising 
in  view  of  the  results  from  Miller  and  Robinson  (1984)  and 
McWilliams  and  Shen  (1960).  A  future  possibility  that  bears  in¬ 
vestigation  is  that  of  embedding  limited-area  models  with  better 
physics,  higher  horizontal  or  vertical  resolution  in  simpler  eddy- 
resolving  basin- scale  models,  or  both.  This  approach  has  en¬ 
joyed  significant  success  in  short-range  atmospheric  forecasts 
(Miyakoda  and  Rosati,  1977),  but  in  situations  where  data  pro- 
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vide  al  least  marginal  resofetkm  of  the  primary  features  at  afi 
levels  in  the  vertical. 

The  executive  summery  of  the  U.S.  Navy-sponsored  Ocean 
Prediction  Workshop  (Mooers  et  a).,  1M2)  slates  that  “the 
needed  capability  of  highest  priority  is  a  portable,  fidl  water  col¬ 
umn,  limited  area,  several-day  to  several-week  prediction  sys¬ 
tem.**  In  the  remainder  of  lids  paper  we  will  discuss  an  alternate 
approach  that  was  not  considered  in  the  Proceedings  of  the 
Ocean  Prediction  Workshop,  but  is  capable  of  overcomint  the 
serious  problems  of  open  boundaries,  computing  power,  and 
data  input  within  the  next  ten  years,  and  meeting  the  description 
of  a  comprehensive  ocean  model  treating  Class  2,  Class  3,  and 
certain  Class  1  phenomena  (mixed  layers  and  equatorial  and 
some  coastal  processes  with  short  time  scales)  listed  in  Table  1. 
The  emphasis  will  be  on  the  treatment  of  Class  2  phenomena, 
such  as  meandering  currents  and  eddies,  because  they  are  the 
most  problematic  in  terms  of  horizontal  resolution,  data  acqui¬ 
sition,  and  data  assimilation,  and  because  satellite  altimeter  data 
should  have  the  greatest  impact  on  this  class.  In  designing  the 
predictive  system,  only  elements  that  are  dearly  necessary  and 
feasible  within  the  next  decade  will  be  considered.  This  does  not 
exclude  the  potential  value  of  other  elements.  Of  necessity  some 
ideas  are  included  that  have  not  been  adequately  tested  and  re¬ 
quire  additional  research. 


A  Cost-Effective  Design  for  Numerical  Ocean  Prediction 
Models 


There  are  two  basic  approaches  that  have  been  used  in  ocean 
model  design .  We  will  call  these  the  ocean  physics  approach  and 
the  ocean  dynamics  approach,  approaches  that  have  counter¬ 
parts  applicable  in  almost  any  field  of  science.  The  ocean  physics 
approach  (1)  concentrates  on  ocean  physics  more  than  ocean 
dynamics  in  designing  the  model,  (2)  makes  the  physics  and  ge¬ 
ometry  of  the  model  as  realistic  as  possible,  and  (3)  assumes  that 
with  good  physics  the  appropriate  ocean  dynamics  will  be  pres¬ 
ent  automatically,  including  some  that  might  not  have  been  an- 
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ticipated.  The  ocean  dynamics  approach  (I)  concentrates  on 
ocean  dynamics  more  than  ocean  physics  m  designing  the  model. 
(2)  designs  a  model  as  simple  and  efficient  (cost-effective)  as 
possible  while  retaining  the  ability  to  represent  the  desired  ocean 
phenomena,  and  (3)  relies  on  investigation  and  understanding  of 
ocean  dynamics  for  successful  application.  Both  the  ocean  phys¬ 
ics  and  die  ocean  dynamics  approaches  require  validation  studies 
to  ensure  that  the  desired  phenomena  are  adequately  repre¬ 
sented. 

As  will  be  shown  in  the  section  on  computing  requirements, 
the  differences  in  the  computational  resources  needed  for  ap¬ 
plication  of  these  two  approaches  are  dramatic,  enough  so  to 
demand  the  use  of  the  ocean  dynamics  approach  for  at  least  the 
next  decade.  Even  when  they  are  computationally  tractable,  the 
complicated  models  resulting  from  the  ocean  physics  approach 
tend  to  be  cumbersome  to  analyze  mid  understand  because  of 
the  large  computational  resources  required  for  each  experiment 
and  because  of  the  huge  number  of  factors  affecting  the  results. 
Furthermore,  as  die  complication  of  the  model  is  increased, 
along  with  the  potential  for  greater  realism,  there  can  be  in¬ 
creased  potential  for  pathological  behavior,  particularly  when  a 
number  of  parameterized  processes  are  included.  Thus  it  is  log¬ 
ical  to  start  with  the  ocean  dynamics  approach  and  build  toward 
more  complicated  models  as  resources  and  understanding  per¬ 
mit,  starting  with  more  realistic  models  of  semienclosed  seas  or 
limited-area  models  embedded  in  simpler  basin-scale  models. 
Even  these  applications  of  an  eddy-resolving  model  derived  from 
the  ocean  physics  approach  are  a  strain  on  current  computer 
capability,  as  can  be  seen  from  the  later  discussion  of  computing 
requirements. 

In  die  following  paragraphs,  the  ocean  dynamics  approach  is 
used  to  design  a  cost-effective  numerical  ocean  prediction  sys¬ 
tem  for  naval  application.  The  thermal  structure  is  the  output 
from  ocean  prediction  models  that  is  currently  of  greatest  in¬ 
terest  to  the  Navy.  Its  mesoscale  variability  is  determined  pri¬ 
marily  by  turbulent  mixing  and  by  the  horizontal  and  vertical 
flow  fields  associated  with  mqjor  current  systems,  eddies, 
oceanic  fronts,  and,  particularly  near  the  equator  and  coastal 
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boundaries,  by  baroctinfc  waves  tad  wind-driven  upwetting  and 
downwelhng.  High  vertical  resolution  (a  few  meters  ia  tbe  upper 
ocean)  is  desired  for  mixed-foyer  modeling  and  for  the  thermal 
field  output  for  naval  operations.  In  contrast,  surprisingly  good 
simulations  of  nugor  current  systems  and  eddies  have  been  ob¬ 
tained  by  using  as  little  as  a  single  vertical  mode  (e.g.,  Hurlburt 
and  Thompson,  1900, 1982).  High  horizontal  resolution  ( ~  10  km) 
is  especially  important  for  these  simulations. 

Great  savings  can  thus  be  obtained  by  separating  the  turbulent 
thermodynamic  problem  from  the  circulation  problem  as  much 
as  possible,  keeping  in  mind  how  the  two  interact.  The  changes 
in  the  thermal  field  resulting  from  mixing  are  communicated  to 
the  circulation  field  primarily  through  the  thermal  wind  relation. 
The  circulation  field  is  communicated  to  the  thermal  structure 
through  the  horizontal  and  vertical  flow  fields. 

Consequently,  we  have  adopted  a  strategy  for  developing  cost- 
effective  ocean  prediction  models  for  the  Navy  based  on  partial 
separation  of  the  hydrodynamic  and  thermodynamic  processes. 
Dynamic  ocean  circulation  models  have  been  developed  that  in¬ 
clude  only  a  small  number  of  foyers  and  incorporate  fairly  crude 
mixed-foyer  physics.  In  general,  two  active  foyers  is  the  mini¬ 
mum  to  permit  coexistence  of  the  pycnodme  and  topography 
and  to  afiow  the  possibility  of  barocimic  instability.  Also,  most 
of  die  energy  in  the  ocean  is  found  m  the  two  lowest  vertical 
modes,  the  barotropic  and  tbe  first  barochnk  (Pochapsky,  1976; 
McWahams,  1976;  Richman  et  al.,  1977;  Ftieri,  1978;  Bernstein 
and  White,  1974).  Schmitz  and  Holland  (1982)  suggest  that  three 
or  more  foyers  may  be  necessary  to  model  the  North  Atlantic. 
For  operational  use  we  suggest  a  minimum  of  three  foyers  with 
linear  stratification  hi  the  central  foyer,  tan  this  is  a  topic  for 
additional  research,  as  is  the  flirtimum  nsnrt  mutation  of  the  ver- 
tied  structure.  The  horizontal  resolution  should  be  sufficient  to 
resolve  muior  currents,  eddies  tooonranhic  finafnm  These 
circulation  models  dura  supply  vertical  and  horizontal  advaction 
information  to  an  N  x  1-D  modal  dim  coasists  of  multiple  (AO 
applications  oft  out  dimensional  (1-D)  model  on  a  suitable  grid. 
The  N  x  l-D  model  has  high  vertical  rtsolutiou  and  a  mote 
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rroftftff  Bfoctim  hMiftaat  l»  farvcMtiat  the  acatHc  Mixed 
layer  in  the  upper  IPO  m.  He  AT  x  14)  Mill  would  provide 
the  thermal  structure  forecast  with  high  vertical  resolution.  To 
iuterfoce  with  the  N  x  14)  Model,  the  layer  depth  everaped  ve¬ 
locities  froM  foe  circulation  model  mm  he  angnmuted  hy  a  per* 
turbatioe  chodatioa  pattsmed  after  Thompson  (1974).  la  the 
peituihatioa  circtdatioa  foe  velocities  aad  messes  ere  Matched 
at  each  layer  boundary  to  provide  dynamically  coasisteat  stooth 
vertical  profiles  of  foe  three  velocity  components,  without  in- 
terpoiatioe.  This  ealetdatioa  iacludfs  foe  effects  of  horizontal 
density  gradients,  vertical  friction  aad  rotation  to  account  for  I 

dynamic  features  such  as  gravitational  circulation,  aad  Ekatan 
drift,  pumping,  and  suction. 

Another  advantage  of  the  partial  decoupling  of  the  turbulent 
thermodynamic  and  the  hydrodynamic  predictions  is  the  one- 
dimensional  nature  of  foe  N  x  1-D  model.  Since  each  vertical 
profile  is  predicted  independently,  the  horizontal  distribution  of 
the  grid  points  is  much  more  flexible  in  both  coven^e  and  spac¬ 
ing.  The  horizontal  resolution  of  foe  N  x  l-D  model  should  be 
sufficient  to  represent  the  forecast  phenomena,  but  it  is  not  dic¬ 
tated  by  numerics  or  dynamics.  Thus,  for  a  given  feature,  the 
grid  spacing  c&  ri  be  two  to  four  times  greater  than  in  a  circulation 
model.  Furthermore,  different  sources  for  flow  fields  can  be  used 
by  different  regions  of  dm  same  N  x  l-D  model.  An  N  x  l-D 
model  from  NORDA  designated  TOW  (Thermodynamic  Ocean 
Prediction  System)  is  already  in  operational  ate  at  the  Fleet  Nu- 
merical  Oceanography  Canter  (Clancy  and  Martin.  1911;  Clancy 
and  Pottak.  1993).  Circulation  rwfrh  are  under  devakxmieat  but 
have  not  yet  been  tested  operationally.  So  for,  they  have  been 
evaraame  pranaruy  oy  mvesnganng  sunutauon  stun,  ocean  oy- 
nantic  properties,  and  sensitivity  to  variations  in  nwdel  param- 
mars,  pqraci,  psumeiry,  lesom&ou,  icrcmg  mnctxms,  ana  m- 
Wal  data.  Studies  of  predictability  and  predictive  akiH  are  also 
necessary.  As  noted  earlier,  seme  work  on  these  topics  has  been 
done  for  open  ecsaa  limited  area  models  (hftSer  and  Robinson, 
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Atmospheric  forecast  Models  are  initialized  from  observations 
that  have  been  analysed  and  dynamically  balanced  at  all  vertical 
levels.  Often,  a  prcvioas  forecast  is  the  first  guess  for  the  analysis 
used  as  foe  initial  state.  The  principal  object  is  to  forecast  waves 
and  eddies  with  scales  of  a  few  thousand  kilometers.  These 
scales  are  marginally  to  adequately  resolved  by  the  observations. 
The  initial  state  is  the  primary  input  data  upon  which  the  forecast 
normally  depends,  except  near  the  surface. 

With  a  few  exceptions  (Haney,  1990;  Clency  and  Martin,  1981; 
Carter  and  Robinson,  1981;  Ocean  Tomography  Group,  1982), 
foe  approach  used  in  meteorology  of  initializing  a  model  from 
observed  data  at  all  vertical  levels  is  not  feasible  in  oceanography 


and  eddies  are  much  smaMer,  8(100  km). 

The  Fleet  Numerical  Oceanography  Center  receives  only 
about  two  hundred  XBT  reports  a  day  for  the  upper  kilometer 
of  foe  global  ocean  (Iodic  and  Fetit,  1982).  However,  seven 
thousand  such  reports  a  day  are  required  to  cover  foe  global 
ocean  with  a  4Mm  grid  once  every  thirty  days  (roughly  equiv¬ 
alent  to  every  three  days  in  the  atmosphere).  MrWjMiami  (1982) 
suagests  294m  resolution  for  moeossals  features.  White  and 
Bernstein  (1979)  suggest  90  km  for  mesoscnlc  eddies  in  the  west¬ 
ern  Fadftc,  but  morale  large  errors.  In  any  case,  whether  the 
iastramsatiti nu  h  XBTs,  air  foropped  XBTs  (AXBTs).  or  sat- 

sivcrofto  witli  foe 

foe  grudlctiou  of  mesoscali  ocean  features,  based  on  con  van 
demdsubsurfoos  data,  is  flaaadafty,  bgMeafty,  and  poiMcaMy 


few  days  (not  every  thirty)  if  we  art  to  perform  even 
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and  short-range  forecasts  of  a  few  days  (Miller  and  Robinson, 
1984). 

Acoustic  tomography  is  a  promising  new  approach  to  sub¬ 
surface  ocean  measurement  that  may  permit  usefbl  coverage  and 
resolution,  because  the  number  of  vertical  profiles  increases 
quadraticaHy  rather  than  linearly  with  the  number  of  deploy¬ 
ments  (Munk  and  Wunsch,  1979,  !962a,b;  The  Ocean  Tomog¬ 
raphy  Group,  1982).  However,  Munk  and  Wunsch  (1982a)  do 
not  suggest  that  eddy-resolving  arrays  will  be  commonplace  or 
that  they  will  cover  a  substantial  portion  of  the  global  ocean  by 
the  1990s.  They  do  propose  some  non-eddy-resolving  arrays 
covering  a  2500-km  scale.  They  also  note  the  important  con¬ 
straint  that  altimeter  data  can  provide  for  the  tomographic  in¬ 
version  problem. 

As  discussed  below,  the  subsurface  data  now  routinely  avail¬ 
able  are  insufficient  to  define  initial  states  that  would  provide 
useful  forecasts  of  mesoscale  features,  such  as  meandering  cur¬ 
rents  and  eddies,  regardless  of  the  methodology  used  for  data 
analysis  and  assimilation.  Assimilation  of  extremely  sparse  sub¬ 
surface  data  by  an  eddy-resolving  forecast  model  could  even  do 
more  harm  than  good.  The  oceanic  subsurface  data  are  so  sparse 
that  in  many  cases  the  deviations  of  the  data  points  from  the 
initial  guess  could  simply  define  spurious  eddies  with  scales 
determined  by  the  influence  fimction  of  the  analysis  scheme, 
eddies  that  could  persist  in  a  several-month  forecast.  For  ex¬ 
ample,  sparse  observations  of  the  Gulf  Stream  could  result  in 
large-amplitude  eddies  in  the  analysis  because  the  observations 
were  insufficient  to  define  a  coherent  and  narrow  meandering 
current. 

Another  problem  with  an  inadequate  initial  state  is  that  some 
of  the  forecast  changes  would  represent  model  spin-up  of  dy¬ 
namic  features  compatible  with  the  physics  and  geometry  of  the 
model,  not  evolution  of  die  ocean  circulation.  For  example,  if  a 
forecast  mode)  resolved  the  Gulf  Stream,  but  the  initial  state  did 
not,  the  model  would  probably  predict  a  narrowing  and  accel¬ 
eration  of  the  Gulf  Stream,  an  event  that  did  not  occur  in  the 
ocean.  It  is  essential  that  both  the  model  and  the  initial  state 


38 


Harley  E.  Hurlburt 


adequately  resolve  the  phenomena  of  interest  in  the  forecast. 
However,  the  circulation  model  needs  horizontal  resolution  two 
to  four  times  finer  than  required  for  adequate  representation  by 
the  data,  in  order  (1)  to  allow  accurate  numerical  integration 
(Grammeftvedt,  1969;  Haidvogel  et  al.,  1980;  Robinson  and 
Haidvogel,  1980)  and  (2)  to  minimize  artificial  damping  of  im¬ 
portant  scales  of  motion  and  associated  nonlinear  interactions 
by  scale-selective  parameterizations  of  subgrid-scale  horizontal 
mixing  (Holland,  1978).  Even  with  adequate  data,  if  the  initial 
state,  the  atmospheric  forcing,  and  the  model  are  independent 
factors  in  the  forecast,  it  is  inevitable  that  model  spin-up  will 
make  some  contribution  to  the  changes  forecast  by  the  model, 
a  phenomenon  that  also  occurs  in  numerical  weather  prediction 
and  is  known  as  the  drift  toward  model  climatology.  This  drift 
can  result  from  both  inadequate  specification  of  the  initial  state 
and  from  differences  between  real  and  model  climatology. 

Even  though  subsurface  data  are  generally  sparse,  where  they 
are  available  with  sufficient  coverage  and  density  to  define  me- 
soscale  features  (nowhere  at  present),  they  should  be  used  in 
eddy-resolving  forecast  models,  particularly  if  they  are  available 
on  a  regular  basis  and  are  sufficient  to  define  space-time  cor¬ 
relation  functions  (see  Bretherton  et  al.,  1976;  Carter  and  Rob¬ 
inson,  1981).  Furthermore,  subsurface  data  must  play  a  crucial 
role  in  the  local  and  regional  validation  of  forecast  models. 

Comeaniam  ef  Surface  km  Subsurface  Data  by  Ocean  Circulation 
Models  and  the  Unigue  Rak  ef  Altimeter  Data 

The  preceding  subsection  has  indicated  that  if  forecasts  of  in¬ 
dividual  current  meanders  and  eddies  are  forced  to  rely  on  in 
situ  data  alone,  then  the  forecasts  will  at  best  be  limited  in  cov¬ 
erage  and  restricted  to  a  short  range,  because  of  the  limited  re¬ 
gions  that  can  be  adequately  sampled.  This  subsection  considers 
die  potential  for  ocean  forecasting,  based  on  observations  from 
satellites  that  have  nearly  global  coverage,  but  only  at  the  sur¬ 
face.  Satellites  and  atmospheric  models  provide  the  only  clear 
prospects  for  fields  with  global  coverage  that  me  useful  opera¬ 
tionally  as  initial  data  or  forcing  functions  for  ocean  forecasts  of 
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meandering  currents,  eddies,  and  frontal  locations.  Satellites 
provide  good  prospects  for  well-observed  fields  of  sea-surface 
elevation,  surface  wind  stress,  sea-surface  temperature,  and  sur¬ 
face  heat  flux  that  could  be  used  in  ocean  prediction. 

With  the  atmospheric  models  and  the  satellite  data,  there  is 
good  potential  for  well-observed  fields  at  the  surface,  which  are 
useful  as  input  to  ocean  circulation  models,  but  these  do  not 
provide  subsurface  data.  At  best,  subsurface  fields  are  ade¬ 
quately  known  only  in  the  mean  plus  a  mean  annual  variation, 
except  possibly  in  limited  areas.  If  a  forecast  model  is  to  rely 
on  data  from  satellites  and  atmospheric  models,  then  well-ob¬ 
served  surface  information  must  be  converted  into  information 
about  the  internal  structure  of  the  ocean,  a  burden  that  might 
logically  be  placed  on  mixed-layer  models  for  the  mixed  layer 
and  on  ocean  circulation  models  for  subsurface  information 
about  current  systems,  eddies,  and  fronts.  So  far  this  is  an  un¬ 
proven  capability  and  an  important  topic  for  research.  However, 
there  are  some  good  hints  that  this  is  feasible. 

One  important  hint  is  the  small  number  of  modes  required  to 
represent  the  vertical  structure  of  the  ocean  (Pochapsky,  1976; 
McWilliams,  1976;  Richman  et  al.,  1977;  Flier],  1978),  a  point 
corroborated  by  Cheney  (1982),  who  has  shown  the  close  rela¬ 
tionship  between  the  surface  dynamic  height  and  the  depth  of 
the  15°C  isotherm  in  the  Sargasso  Sea.  This  point  is  illustrated 
by  a  two-layer  reduced-gravity  model  that  consists  of  a  single 
internal  vertical  mode,  the  simplest  ocean  circulation  model  with 
potential  usefulness  in  simulation  and  prediction.  The  model  con¬ 
tains  an  active  upper  layer  and  a  lower  layer  that  is  infinitely 
deep  and  at  rest.  The  interface  between  the  layers  represents 
the  pycnocline,  the  sharp  change  in  density  with  depth  that  oc¬ 
curs  over  most  of  the  world  ocean. 

One  notable  feature  of  this  model  is  its  capability  for  initial¬ 
ization  by  altimeter  data  alone,  except  near  the  equator  where 
knowledge  of  the  wind  stress  is  important  for  initialization  in 
addition  to  forcing.  A  second  feature  is  its  one-to-one  corre¬ 
spondence  between  the  variations  in  the  sea-surface  elevation 
related  to  ocean  currents  and  the  variations  in  the  depth  of  the 
pycnocline,  and  a  third  feature  is  its  remarkable  ability  to  sim- 
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ulate  certain  ocean  phenomena.  This  ability  has  been  demon¬ 
strated  for  a  variety  of  phenomena  in  several  different  regions, 
for  example,  the  equatorial  Pacific  (Kindle,  1979;  Busalacchi  and 
O'Brien,  1980;  Busalacchi  et  al.,  1983),  the  Indian  Ocean  (Cane, 
1980),  the  Somali  Current  (Lin  and  Hurlburt,  1981),  the  Gulf  of 
Mexico  (Hurlburt  and  Thompson,  1980,  1982),  ocean  response 
to  a  hurricane  (O'Brien  and  Reid,  1967;  Chang  and  Anthes,  1978; 
Price,  1981),  the  Gulf  of  Guinea  upwelling  (Adamec  and  O'Brien, 
1978),  the  Alboran  Sea  (Preller  and  Hurlburt,  1982)  and  the  west¬ 
ern  Mediterranean  Sea  (Preller  and  Heburn,  1983).  (A  few  of 
these  models  have  included  more  than  one  internal  vertical  mode 
or  mixed-layer  physics.)  The  preceding  references  and  the  spe¬ 
cial  compatibility  with  altimeter  data  offer  hints  that  a  reduced- 
gravity  model  with  a  single  vertical  mode  may  be  able  to  forecast 
a  substantial  part  of  the  potentially  predictable  oceanic  varia¬ 
bility  associated  with  major  current  systems  and  mesoscale  ed¬ 
dies. 

However,  with  only  a  single  (internal)  vertical  mode,  reduced 
gravity  models  cannot  be  recommended  for  genera)  use  as  op¬ 
erational  eddy-resolving  forecast  models,  because  they  do  not 
account  for  the  potentially  important  effects  of  topography  and 
baroclinic  instability.  The  simplest  ocean  model  that  can  include 
these  effects  is  a  two-active-layer  model  with  an  internal  (bar- 
oclinic)  mode  and  an  external  (barotropic)  mode,  the  latter  rep¬ 
resenting  the  depth-independent  flow.  The  importance  of  having 
at  least  two  vertical  modes  for  the  Gulf  Stream  has  been  dem¬ 
onstrated  by  the  work  of  Holland  and  Lin  (1975),  Rhines  (1977), 
Semtner  and  Holland  (1978),  and  Hurlburt  and  Thompson  (1984), 
and  for  the  Antarctic  Circumpolar  Current  by  McWilliams  et  al. 
(1978).  Schmitz  and  Holland  (1982)  suggest  that  three  or  more 
layers  may  be  required  for  the  North  Atlantic.  It  is  also  note¬ 
worthy  that  the  sea-surface  elevation  in  layered  models  is  bal¬ 
anced  by  the  depth-averaged  geostrophic  velocity  in  the  upper¬ 
most  layer,  while  the  oceanic  current-related  sea-surface 
elevation  is  balanced  by  the  geostrophic  component  of  surface 
currents.  Thus  the  models  tend  to  underestimate  the  amplitude 
of  the  current-related  sea-surface  elevation,  a  situation  that 
should  be  taken  into  account  when  one  is  using  altimeter  data 
in  ocean  models. 
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Munk  and  Wunsch  (1982a)  corroborate  the  severe  limitations 
of  conventional  in  situ  data  (and  the  promise  of  acoustic  to¬ 
mography)  for  describing  the  world  ocean.  In  addition,  they 
cite  the  importance  of  satellite  data  at  the  surface  and  the  im¬ 
portance  of  models  in  converting  these  surface  data  into  sub¬ 
surface  information. 


Simulation  Skill  of  a  Low-Vertical  Resolution  Model 
for  the  Gulf  of  Mexico 

As  an  example.  Figure  2  illustrates  the  simulation  skill  of  a  model 
for  the  Gulf  of  Mexico  with  two  active  layers,  realistic  coastline 
geometry  and  bottom  topography,  and  horizontal  resolution  suf¬ 
ficient  for  major  current  systems  and  eddies.  Surprisingly,  no 
eddy-resolving  ocean  modeling  study  published  until  now  has 
included  all  of  these  features  in  a  deep  ocean  basin  where  plan¬ 
etary  vorticity  advection  was  important  and  the  model  was  in¬ 
tegrated  to  statistical  equilibrium.  This  calculation  was  per¬ 
formed  by  A.  Wallcraft  of  JAYCOR,  Alexandria,  VA  (personal 
communication,  1982)  using  the  model  of  Hurlburt  and  Thomp¬ 
son  (1980).  The  model's  ability  to  handle  detailed  coastline  ge¬ 
ometry  was  Wallcraft 's  addition.  The  model  was  driven  from 
rest  to  statistical  equilibrium  solely  by  a  steady  inflow  through 
the  Yucatan  Straits  between  Cuba  and  Mexico,  which  was  com¬ 
pensated  by  outflow  through  the  Florida  Straits.  The  model  pa¬ 
rameters  are: 

Upper/lower  layer  inflow  transport  *  26/4  m3/s 
Horizontal  eddy  viscosity  A  =  300  m2/s 
The  Coriolis  parameter  at  the  southern  boundary  /„  =  5  x 
10"5  s-1 

Gravitational  acceleration  g  =  980  cm/s2 
Reduced  gravity  g'  =  3(tf,  +  H2)/H2  cm/s2 
Reference  thicknesses  of  the  layers  Ht  -  200  m,  H2  -  3400 
m,  p  *  dfldy  *  2  x  10-13  cm' ’ s" 1 
Surface  and  interfacia]  stresses  =  0 
Stress  coefficient  for  quadratic  bottom  stress  C*  =  2  x  10“ 3 
Horizontal  grid  increments  for  each  variable  Lx  -  Ay  =  25 
km 


FIGURE  2.  (a)  Instantaneous  view  of  the  interface  deviation  in  a  two-layer 
simulation  of  the  Gulf  of  Mexico  driven  from  rest  to  statistical  equilibrium 
solely  by  inflow  through  the  Yucatan  Straits.  The  contour  interval  is  25  m, 
with  solid  contours  representing  downward  deviations.  ( b )  Depth  of  the  22°C 
isothermal  surface.  4-18  August  1966  ( Alaminos  cruise  66-A-l  I),  from  Leipper 
(1970).  The  contour  interval  is  25  m. 
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Time  step  Ar  *  1  h 

Inflow  spin-up  time  constant  *  thirty  days 

Minimum  depth  of  the  bottom  topography  *  500  m 

Figure  2  compares  “instantaneous"  upper  ocean  flow  patterns 
(a)  from  the  numerical  model  and  (b)  from  observations  by  Leip- 
per  (1970).  The  Loop  Current  is  the  mqor  current  system 
depicted.  At  this  point  it  is  about  to  shed  an  eddy.  The  Loop 
Current  is  observed  to  penetrate  the  Gulf  of  Mexico,  bend 
westward,  and  shed  an  eddy  with  a  period  of  approximately  one 
year.  This  “annual"  cycle  was  long  thought  to  be  due  tc  seasonal 
variations  in  the  flow  through  the  Yucatan  Straits  (Cochrane, 
1965).  However,  the  model  Loop  Current  exhibits  an  approxi¬ 
mately  annual  eddy-shedding  period  when  the  inflow  is  steady, 
contradicting  the  earlier  hypothesis.  Although  time  variations 
are  not  essential,  they  can  still  play  a  significant  role  in  the  eddy 
shedding  (Huriburt  and  Thompson,  1980).  The  model  Loop  Cur¬ 
rent  also  spontaneously  shed  eddies  with  realistic  diameters,  am¬ 
plitudes,  and  westward  propagation  speeds.  The  dynamic  basis 
for  the  agreement  between  the  observations  and  a  circulation 
model  with  low  vertical  resolution  is  discussed  by  Huriburt  and 
Thompson  (1980,  1982). 

After  an  eddy  is  shed  by  the  Loop  Current,  it  propagates  west¬ 
ward  (Figure  3a).  In  this  case  there  was  spontaneous  develop¬ 
ment  of  a  counter-rotating  vortex  pair  when  the  eddy  reached 
the  western  gulf  (Figure  3 b),  a  structure  repeatedly  observed  in 
that  region  (Figure  4).  The  roles  of  the  wind  and  the  eddies  from 
the  Loop  Current  in  the  formation  of  this  structure  have  been  a 
matter  of  some  controversy  (Merrell  and  Morrison,  1981).  Wind 
forcing  was  omitted  from  the  numerical  experiment  shown  in 
Figures  2  and  3.  Although  wind  forcing  is  not  essential,  a  sig¬ 
nificant  role  for  it  has  not  been  ruled  out. 

Strategy  for  the  Assimilation  of  Surface  Data  in 
Ocean  Circulation  Models 

We  have  already  noted  the  simulation  skill  of  models  with  low 
vertical  resolution  and  the  one-to-one  correspondence  between 
the  sea-surface  elevation  and  the  internal  structure  in  a  model 
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FIGURE  3.  (a)  Interface  deviation  from  the  Gulf  of  Mexico  simulation  at  day 
1970  alter  an  eddy  has  separated  from  the  model  Loop  Current  and  propagated 
westward,  (b)  Ninety  days  later  the  mqor  anticydonic  eddy  at  day  1970  has 
developed  into  a  counter-rotating  vortex  pair  in  the  western  gulf.  The  cyclonic 
vortex  is  to  the  north  and  the  anticylonic  to  the  south.  The  contour  interval 
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FIGURE  4.  Counter-rotating  vortex  pair  in  the  western  Gulf  of  Mexico  as 
shown  by  the  depth  of  the  I5*C  isotherm  (in  meters),  observed  in  April  1978. 
The  cyclonic  vortex  is  to  the  north  and  the  anticylonic  to  the  south.  The  con¬ 
tour  interval  is  25  m  (from  Merrell  and  Morrison.  1981). 
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that  represents  the  single  vertical  node  explaining  the  majority 
of  the  oceaafc  variance.  In  view  of  this  and  the  discussion  earlier 
in  this  sectiei,  it  seems  preferable  to  rely  on  the  simulation  skill 
of  a  more  realistic  forecast  model  to  distribute  the  potentially 
abundant  surface  information  among  the  small  number  of  vertical 
modes  than  tk>  assimilate  subsurface  in  situ  data  m  regions  where 
they  are  too  aparse  to  adequately  define  mesoscale  features. 

Thus  the  following  procedure  is  recommended:  The  ocean  cir¬ 
culation  model  for  Class  2  and  Gass  3  phenomena  (Table  1) 
should  'be  spun  up  to  statistical  equilibrium  by  using  the  best 
historical  data  available  as  forcing  functions.  This  would  provide 
the  initial  state  for  ocean  forecasting.  The  initial  state  could  be 
updated  by  integrating  the  model  forward  in  time  as  new  data 
for  the  forcing  functions  arrived.  Once  they  became  available, 
scatterometer-derived  wind  stress  data  and  the  current-related 
component  of  the  sea-surface  elevation  from  altimeter  data  could 
be  introduced  as  forcing  functions.  This  should  greatly  improve 
the  accuracy  of  the  initial  state  for  the  ocean  forecast,  and  for 
the  first  time*make  it  possible  to  predict  the  evolution  of  indi¬ 
vidual  current  meanders  and  eddies  on  a  global  basis.  Introduc¬ 
tion  of  the  altimeter  data  would  require  special  care  to  avoid 
serious  unphysical  dynamic  imbalances  in  the  model.  This  ap¬ 
proach,  though  attractive,  is  unproven  at  this  time  and  is  an 
important  tope  for  research.  In  regions  where  adequate  sub¬ 
surface  data  are  available,  they  should,  of  course,  be  used. 
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The  horizontal  and  temporal  resolution  required  for  the  forcing 
functions  is  dictated  primarily  by  the  space  and  time  scales  of 
the  atmosphere.  Thus,  scatterometer-derived  wind  stress  data 
would  be  needed  on  a  daily  basis  to  resolve  the  evolution  of 
atmospheric  storm  systems.  The  spatial  resolution  of  the  wind 
stress  data  should  be  sufficient  (—50  km)  to  resolve  the  wind 
stress  curl  associated  with  atmospheric  cyclones,  anticyclones, 
and  fronts.  Winds  as  close  to  coastlines  as  possible  are  highly 
desirable.  In  these  regions  horizontal  resolution  should  be  20  km 
or  better.  Since  a  satellite  scatterometer  with  an  orbit  altitude 
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of  about  830  km  has  about  a  15004m  swath  width  (with  gaps), 
a  single  satellite  in  coqjunction  with  atmospheric  models  and  the 
existing  data  base  could  meet  these  requirements  adequately 
(Satellite  Surface  Stress  Working  Group.  1982). 

Cheney  and  Marsh  (1981)  have  already  demonstrated  the  abil¬ 
ity  of  the  SEASAT  altimeter  to  detect  the  Gulf  Stream  and  me- 
soscale  eddies  when  using  either  repeat  tracks  or  a  geoid,  and 
Cheney  et  al.  (1983)  have  used  nine  sets  of  three-day  repeat 
tracks  from  SEASAT  to  produce  a  global  map  depicting  the  root- 
mean-square  (RMS)  mesoscale  variability  of  the  sea-surface  el¬ 
evation  during  September  and  October  1978.  Tapley  et  al.  (1982) 
and  Wunsch  and  Gaposchkin  (1980)  review  the  problem  of  ob¬ 
taining  accurate  measurements  of  current-related  variations  in 
the  sea-surface  elevation  from  a  satellite  altimeter  and  the  suc¬ 
cess  achieved  with  SEASAT.  To  be  usdful  in  ocean  circulation 
prediction,  the  altimeter  data  must,  of  course,  resolve  mqjor  cur¬ 
rent  systems  and  eddies  in  time,  amplitude,  and  horizontal  di¬ 
mension.  However,  the  amplitude  dependence  on  horizontal 
scale  should  also  be  noted.  For  currents  in  geostrophic  balance, 

(  x  fv,  - 

where  /  is  the  Coriolis  parameter,  v,  is  the  velocity  of  the  geo¬ 
strophic  surface  current,  g  is  the  acceleration  due  to  gravity,  and 
t)  is  the  sea-surface  elevation  above  the  geoid  related  to  oceanic 
surface  currents.  TLie  maximum  velocity  V  for  typical,  signifi¬ 
cant,  persistent,  mesoscale  (~50  to  500  km)  features  in  the  ocean 
ranges  from  about  10  cm/s  to  1  m/s.  Up  to  2  m/s  is  not  uncommon 
in  the  stronger  currents.  For  an  altimeter  that  can  useftilly  mea¬ 
sure  oceanographic  features  with  a  minimum  amplitude  of  n  » 
10  cm.  the  minimum  radius  of  an  eddy  or  half-width  of  a  current 
L  with  V  -  10  cm/s  is  L  -  100  km  at  AT  latitude  and  L  -  ISO 
km  at  27*  latitude.  For  V  -  100  cm/s,  L  «  10  km  at  42*  latitude, 
and  die  Rossby  number  R„  «  VlfL  -  1.  For  R.  <  1  the  flow  is 
geostrophkally  balanced.  For  Ra  ►  1  vortices  can  develop  cy- 
clostropliic  balance. 
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where  v%  is  die  tangential  velocity  and  r  it  tlw  radius  of  curvature, 
la  this  can  the  wnplfcade  of  a  associated  with  the  vortex  is 
depsadsonaa  v*  but  is  nearly  iadepeadeat  of  the  scale  of  the 
vortex.  Hdwtver,  a  velocity  scale  greater  thaw  2  m/s  is  required 
for  a  signal  ia  the  sea  sarfoce  elevation  greater  than  10  cm. 

From  this  w«  conclude  that  a  sateKte  altimeter  able  to  measure 
variations «i  foe  sea-surfece  ekvatiea  greater  than  5  cm,  which 
are  associated  with  oceanic  currents,  should  provide  usable  in- 
fonuatioa  for  the  more  significant  currents  and  eddies  with  scales 
greater  than  a  few  tens  of  kilometers.  This  accuracy  is  feasible 
with  current  technology  (TOPEX  Science  Working  Group, 
1981).  Altimeter  track  spacing  of  about  30  km  and  at  most  50  km 
at  midfelitudes  is  needed  for  horizontal  resolution  of  most  de¬ 
tectable  eddies.  Experiments  by  Kindle  et  a).  (1984)  have  shown 
that  three  ascending  or  descendu*  tracks  with  error-free  data 
are  sufficient  to  adequately  map  an  eddy  with  a  nearly  circular 
shape.  Repeat  intervals  up  to  a  month  should  provide  adequate 
temporal  nesohition  for  most  mesoscale  features  (Figure  1).  par¬ 
ticularly  with  the  aid  of  predictive  models  to  fill  in  the  temporal 
gaps.  To  minimize  the  time  required  for  altimeter  tracks  to  span 
a  given  eddy,  the  orbit  should  be  chosen  so  that  adjacent  tracks 
are  one  day  apart.  This  requirement  should  override  any  concern 
about  tidy  abasing,  since  the  tides  are  relatively  large-scale  phe¬ 
nomena  and  because  removal  of  foe  tidal  contributions  to  the 
sea-surfrcc  elevation  appears  feasible  (Schwiderski,  1980;  Cart¬ 
wright  and  Alcock,  1981;  Rtowa  and  Hutchinson,  1981;  Dia¬ 
mante  and  Nee,  1981;  liaziegs,  1983).  However,  as  noted  in 
some  of  foe  preceding  references,  accurate  determination  of  foe 
tidal  aliases  end  amnlitudes  would  be  facilkatad  bv  at  once 
having  an  altimeter  in  an  orbit  relatively  free  of  tidal  aliases  with 
loan  periods. 

A  single  mieBite  in  a  twenty-day  repeat  orbit  that  carried  a 
multibeam  altimeter  (see  Bush  et  aJ.,  1984,  in  this  issue)  could 
easily  maet  foe  preceding  requirements.  Minimal  satisfaction  of 
the  requirements  could  be  achieved  by  using  a  single  satellite 
with  a  nadir  beam  altimeter  in  a  forty-day  repeat  orbit,  if  a  nu¬ 
merical  model  was  able  to  bridge  the  temporal  gaps.  Because 
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crabfc  to  reduce  the  temporal  resohstioa  rather  thee  the  ipatial 
resolutioo.  la  the  ahteace  of  aa  adequate  paid,  it  would  he 
necessary  to  use  a  aeaa  sea  sarfhee.  to  that  oaly  deviations  from 
the  aieaa  would  be  available  to  the  model,  la  that  erne  the  model 
could  obtaia  the  meaa  from  its  owe  or  observed  climatology. 
Bandpass  filtering  could  be  required  to  remove  short  (<0  (10 
km))  and  long  wave  (greater  than  a  few  thousand  kilometers) 
errors,  the  scales  with  the  greatest  errors  (TOPEX  Science 
Working  Group,  1901). 


Once  a  particular  ocean  model  design  has  been  chosen  and  the 
universal  need  for  approximately  10-km  grid  resolution  of  ocean 
currents  aad  eddies  has  been  noted,  then  there  is  a  simple  re¬ 
lationship  between  basin  size  and  the  computer  power  required. 
Based  oa  a  three-layer  version  of  the  circulation  model  discussed 
in  the  section  on  cost-effective  design.  Table  2  lists  oceanic  do¬ 
mains  versus  computer  power  required.  The  Class  7  machines 
expected  on  the  market  in  the  middle  and  late  1980s  will  be  the 
first  that  are  adequate  for  all  nuqor  ocean  basins. 

In  view  of  the  huge  computer  power  required,  possibilities  for 
reducing  it  merit  consideration.  The  limited  utility  of  a  model 
with  low  horizontal  resolution  has  already  been  noted  in  the  sec¬ 
tion  on  time  aad  space  scales  aad  Figure  1.  Three  other  possi¬ 
bilities  are  (1)  alternate  model  designs,  (2)  the  use  of  virtual/ 
auxiliary  arnaMry,  aad  (3)  the  use  of  nested  grids  or  subdomains 
of  major  ocean  basins. 

Table  3  compares  one  model  (NORDA's)  developed  by  using 
the  ocean  dynamics  approach  (see  the  beginning  of  the  section 
oa  cost-effective  design)  with  two  models  closer  in  design  to  die 
ocean  physics  approach  (Science  Applications,  lac.  (SA1)  aad 
Dyaalysis).  The  SA1  (Roberts  et  al.,  1980)  aad  Dynalysts  (Bhun- 
berg  aad  Meflor,  1983)  models  me  both  multilevel  models  with 
a  terrain-following  o-coordinate  system.  Approximately  fifteen 
levels  are  the  mmtmum  required  for  adequate  resolution  of  the 
vertical  structure  mid  the  relatively  sophisticated  parameteri¬ 
zation  of  the  aiixed-layer  physics.  The  NORDA  model  can  work 
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with  much  lower  vertical  resolution  by  using  natural  layers  in 
the  oceanic  structure  and  a  Lagrangtan  vertical  coordinate  in 
which  the  layers  change  thickness  and  move  vertically  in  re¬ 
sponse  to  the  ocean  circulation  and  a  crude  parameterization  of 
vertical  mixing. 

In  the  approach  discussed  in  the  section  on  cost-effective  de¬ 
sign,  a  rirculation  model  with  low  vertical  resolution  fed  hori¬ 
zontal  and  vertical  advectioa  information  to  an  N  x  |-D  mixed- 
layer  model  with  hi0i  vertical  resolutioo  and  more  sophisticated 
mixed  layer  physics,  lecauie  the  calculations  at  the  individual 
grid  points  are  independent  m  the  N  x  l-D  model,  the  numerics 
and  the  dynamics  do  not  dictate  horizontal  grid  resolution  finer 
than  regained  to  represent  the  features  of  interest,  and  the  grid 
spacing  can  he  turn  to  four  times  greater  than  in  the  circulation 
model.  Also,  open  bounderiei  are  not  a  problem  and  the  cov¬ 
erage  and  resolution  are  much  more  flexible.  Furthermore, 
mixed  layers  we  Class  1  phenomena  (Table  I)  which  we  re¬ 
stricted  to  short-rangi  forecasts.  Thus  the  N  x  l-D  model  can 
be  e  relatively  inexpensive  calculation  and  is  not  included  in 
Table  J. 
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The  comparison  in  Table  3  is  based  on  estimated  resources 
required  for  a  two- month  forecast  using  an  eddy-resolving  North 
Atlantic  prediction  model.  Two- month  forecasts  are  appropriate 
for  Class  2  phenomena  such  as  meandering  currents  and  eddies. 
These  three  models  were  partly  chosen  to  facilitate  the  com¬ 
parison,  because  all  three  have  executed  smaller  computations 
on  the  same  vector  computer,  the  two-pipe  Texas  Instruments 
Advance*  Scientific  Computer  (Tl/ASC)  at  the  Naval  Research 
Laboratory  in  Washington,  DC.  The  Tl/ASC  times  estimated  for 
the  North  Atlantic  were  reduced  by  1/7.5  to  provide  estimates 
for  a  two-pipe  CYBER  205,  the  most  powerful  computer  cur¬ 
rently  used  for  ocean  modeling  (Geophysical  Fluid  Dynamics 
Laboratory.  Princeton.  NJ;  NASA/Goddard,  Greenbelt,  MD; 
Fleet  Numerical  Oceanography  Center  (one-pipe),  Monterey, 
CA).  The  nominal  execution  speeds  for  the  NOR  DA  model  are 
20  megaflops  on  the  two-pipe  Tl/ASC,  80  megaflops  on  a  one- 
pipe  CYBER  205,  and  150  megaflops  on  a  two-pipe  CYBER  205 
(A.  Walkraft.  personal  communication,  1983). 

The  mdidels  are  listed  at  two  different  resolutions  to  make  two 
points.  (I)  The  SAI  and  Dynalysis  models  require  higher  hori¬ 
zontal  resolution  to  represent  a  given  topographic  feature  with 
the  same  accuracy  as  the  NORDA  model.  The  topography  ap¬ 
pears  undifferentiated  in  the  equations  of  the  NORDA  model, 
but  is  differentiated  in  the  equations  for  the  other  two;  see,  for 
example,  Rhines's  (1977)  discussion  of  the  effects  of  differen¬ 
tiation  on  spectrum  centroids.  (2)  With  double  the  horizontal 
resolution,  the  NORDA  model  stHI  requires  smaller  resources 
than  the  other  two.  The  relative  merits  of  increasing  horizontal 
and  vertical  resolution  should,  of  course,  be  tested,  but  on  a 
smaller  domain  than  the  North  Atlantic.  This  paper  has  advo¬ 
cated  die  former  for  most  Class  2  and  Class  3  phenomena  listed 
in  Table  1,  but  high  vertical  resolution  is  required  when  sophis¬ 
ticated  mixed-layer  physics  are  included.  In  any  case,  the  SAI 
and  Dynalysis  models  appear  intractable  as  eddy-resolving 
models  of  the  larger  ocean  basins  at  least  for  the  next  decade. 
For  a  given  model,  the  memory  requirements  depend  inversely 
on  the  square  and  the  computer  time  inversely  on  the  cube  of 
the  horizontal  resolution. 
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All  three  models  use  primitive  equations  and  retain  a  free  sur¬ 
face.  The  NORDA  and  SA1  models  increase  the  time  step  per¬ 
mitted  by  using  implicit  numerical  schemes  for  at  least  some  of 
the  terms.  However,  the  high  degree  of  implicitization  in  the  SA1 
model  makes  it  inherently  the  most  complicated  and  most  dif¬ 
ficult  of  the  three  to  modify.  The  Dynalysis  model  uses  an  ex¬ 
plicit  numerical  scheme  in  the  horizontal,  but  separates  the  ex¬ 
ternal  from  the  internal  modes.  This  separation  permits  a  much 
longer  time  step  for  the  internal  modes,  because  the  internal 
gravity  wave  speed  is  much  less  than  the  external.  The  time  step 
limitation  in  the  NORDA  model  differs  from  that  of  a  quasi- 
geostrophic  model  only  by  a  1/|/|  stability  criterion  imposed  by 
the  numerical  scheme  for  the  Coriolis  force.  At  the  resolution 
used  for  the  North  Atlantic  model  in  Table  3,  the  quasigeo- 
strophic  approximation  would  have  little  effect  on  the  computer 
time  required.  Although  this  approximation  can  be  advantageous 
in  some  situations,  it  was  not  used  in  the  NORDA  model  because 
it  is  not  valid  near  the  equator,  for  large-amplitude  topography, 
large  scales  (*>1 ,000  km),  or  large  vertical  excursions  of  the  pyc- 
nocline.  For  example,  Hurlburt  and  Thompson  (1984)  found  that 
the  effects  of  large-  and  small-amplitude  seamounts  on  the  Gulf 
Stream  were  quite  different. 

The  use  of  auxiliary  memory  and  asynchronous  buffering  has 
long  been  a  popular  practice  for  atmospheric  forecast  models 
and  was  used  in  early  Tl/ASC  versions  of  the  Dynalysis  ocean 
model  (Blumberg  and  Mellor,  1981).  (However,  the  times  in 
Table  3  are  based  on  a  core-contained  version  [Blumberg  and 
Mellor,  1984],  as  for  the  other  two  models.)  A.  Wallcraft  has 
noted  (personal  communication,  1982)  that  this  procedure  is  not 
feasible  for  ocean  prediction  models  unless  the  input/ouiput 
speed  in  megawords/s  for  virtual  or  auxiliary  memory  is  0. 1  to 
0.3  times  the  usable  megaflop  rate.  This  assumes  spending  equal 
time  for  model  execution  and  memory  swapping  input/output. 
On  the  CYBER  205  computer,  time-dependent  ocean  models 
must  be  core-contained  because  swapping  data  in  and  out  of 
central  memory  is  too  slow  in  relation  to  the  central  processor 
speed,  by  a  factor  of  about  a  hundred.  Even  though  this  situation 
may  change  (for  example,  by  using  solid  state  secondary  storage 
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devices  or  hierarchical  memory  systems  that  effectively  provide 
buffering  between  disks  and  main  memory),  central  processor 
speeds  of  ~  1  gigaflop  are  still  required  for  the  larger  ocean  ba¬ 
sins. 

Ocean  predictions  on  subdomains  of  major  ocean  basins  are 
the  third  and  superficially  most  attractive  way  to  reduce  the  com¬ 
putational  power  required.  As  discussed  in  earlier  sections,  ex¬ 
tensive  open  boundaries  are  the  primary  difficulty  with  this  ap¬ 
proach.  Open  boundaries  can  be  tolerated  in  situations  where 
most  of  the  interesting  phenomena  are  generated  within  the  do¬ 
main  or  in  cases  where  the  boundary  information  required  is 
fairly  simple,  such  as  flow  through  a  strait  or  total  transport 
through  the  boundary.  Thus,  semienclosed  domains  with  small 
open  boundary  segments  or  major  ocean  basins  with  distant  open 
boundaries  are  tractable. 

In  sub-basin-scale  domains  where  open  boundaries  must 
admit  significant  meandering  current  systems  and  eddies,  the 
open  boundaries  become  the  most  serious  limitation  to  data  as¬ 
similation  and  forecasting  skill,  unless  accurate  surface  and  sub¬ 
surface  boundary  information  on  the  currents  and  eddies  is  avail¬ 
able  for  the  duration  of  the  forecast  (Miller  and  Robinson,  1984). 
The  requirement  for  information  at  the  open  boundaries  becomes 
even  more  severe  when  the  role  that  models  must  play  in  assim¬ 
ilating  data  and  performing  nowcasts  is  considered.  As  noted 
earlier,  the  conversion  of  well-observed  surface  information  into 
subsurface  information  by  means  of  circulation  models  is  not  an 
instantaneous  process,  except  in  models  with  a  single  vertical 
mode.  Thus  even  if  adequate  data  were  available  at  the  surface, 
errors  at  the  subsurface  open  boundaries  would  propagate  into 
the  domain  during  the  data  assimilation  process,  degrading  the 
ability  to  initialize  the  model  and  make  even  short-range  fore¬ 
casts  of  meandering  currents  and  eddies.  This  situation  is  even 
further  aggravated  when  the  open  domain  of  interest  moves  with 
time,  as  expected  in  many  Naval  applications.  In  basin-scale 
domains,  assimilation  of  altimeter  data  during  the  nowcasting 
phase  should  help  to  prevent  the  spread  of  errors  introduced  at 
artificial  boundaries  more  than  a  yet-undetermined  distance  into 
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the  mode!  domain.  A  small  amount  of  crucially  sited  in-situ  data 
should  also  help. 

Given  a  limited-area  model,  three  possible  ways  to  reduce  the 
open  boundary  problem  are  (1)  frequent  (every  few  days)  sub¬ 
surface  sampling  on  the  boundaries  (Robinson  and  Haidvogel. 
1980);  (2)  having  a  satellite  altimeter  in  a  few-day  repeat  orbit, 
using  the  repeat  orbits  as  domain  boundaries,  and  empirically 
converting  the  surface  into  subsurface  information;  and  (3)  ob¬ 
taining  boundary  information  from  an  eddy-resolving  basin-scale 
model.  Data  with  less  temporal  resolution  are  required  for  the 
interior.  The  first  method  is  impractical  for  wide  application 
within  the  next  decade.  The  second  is  untested  and  does  not 
provide  data  for  the  interior  of  the  domain.  The  results  from 
Miller  and  Robinson  (1984)  suggest  that  the  success  of  the  third 
method  would  depend  strongly  on  the  quality  of  the  boundary 
conditions  from  the  basin-scale  model.  In  the  first  two  methods, 
the  accuracy  of  the  boundary  information  at  future  times  is  likely 
to  limit  the  forecasts  to  a  few  days,  although  statistical  forecasts 
of  the  boundary  conditions  may  aid  significantly  in  some  situ¬ 
ations  (Robinson  and  Tu,  1982). 

Forecasts  of  up  to  a  few  months  are  potentially  permitted  by 
the  dynamics  of  phenomena  in  Classes  2  and  3  (Table  1),  phe¬ 
nomena  such  as  meandering  currents  and  eddies.  Forecasts  of 
this  duration  would  be  highly  advantageous  in  some  Naval  ap¬ 
plications,  since  they  would  greatly  facilitate  dissemination  of 
information  to  a  fleet,  permit  advance  planning  when  knowledge 
of  the  ocean  environment  was  important,  and  allow  substantia] 
knowledge  of  the  ocean  circulation  and  thermal  structure  up  to 
a  few  months  after  the  loss  of  satellites  or  other  sources  of  data. 

In  most  real-time  circulation  model  applications  to  open  do¬ 
mains,  models  with  basin-scale  domains  provide  the  most  prac¬ 
tical  source  of  data  for  initialization  and  open  boundary  condi¬ 
tions,  conditions  required  at  all  vertical  levels  and  for  the 
duration  of  the  forecast.  However,  if  the  basin-scale  model  is  to 
provide  useful  information,  it  too  must  resolve  the  meandering 
currents  and  eddies  in  the  forecast.  Thus  we  cannot  escape  the 
need  for  high-resolution  (—10  km)  models  with  domains  covering 
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major  ocean  basins,  the  need  for  satellite  altimeter  and  scat- 
terometer  data  with  global  coverage,  and  the  need  for  Class  7 
computers  to  perform  the  numerical  integrations  for  the  fore¬ 
casts,  for  spin-up  of  the  subsurface  initial  state,  for  model  val¬ 
idation  studies,  and  for  studies  of  model  dynamics  and  their  re¬ 
lation  to  the  ocean.  The  subsurface  initialization,  the  model  val¬ 
idation  studies,  and  the  studies  of  model  dynamics  can  require 
multi-year  integrations  to  allow  the  model  to  reach  statistical 
equilibrium^ 

Summary  and  Conclusions 

Within  the  next  decade  there  is  high  potential  for  a  global  real¬ 
time  description  and  prediction  capability  for  oceanic  eddies, 
meandering  current  systems,  and  frontal  positions  on  time  scales 
up  to  a  few  months  (despite  atmospheric  predictability  of  only 
about  one  week;  see  Table  1).  The  critical  elements  required  to 
realize  this  potential  are  (1)  adequate  data  input,  (2)  adequate 
computing  power,  and  (3)  properly  designed  and  adequately  val¬ 
idated  ocean  models  for  both  data  assimilation  and  prediction. 

Highly  efficient  models  with  low  vertical  resolution  have  dem¬ 
onstrated  a  remarkable  ability  to  simulate  meandering  currents 
and  eddies,  but  a  high  horizontal  resolution  (~  10  km)  is  required 
to  resolve  them.  To  avoid  serious  limitations  on  data  assimilation 
and  predictive  skill,  it  is  advantageous  to  use  domains  that  cover 
major  ocean  basins  or  semienclosed  seas  where  the  phenomena 
of  interest  are  generated  within  the  domain,  and  where  the  open 
boundary  segments  are  either  small  (such  as  flow  through  a 
strait)  or  very  distant  from  the  region  of  interest.  Circulation 
models  on  subdomains  of  m^jor  ocean  basins  that  have  extensive 
open  ocean  boundaries  are  extremely  dependent  on  some  source 
like  a  large-scale  model  for  boundary  conditions  during  both  the 
data  assimilation  and  prediction  phases.  To  provide  useful 
boundary  conditions,  the  large-scale  model  must  also  have  res¬ 
olution  sufficient  for  the  currents  and  eddies.  To  model  m^jor 
ocean  basins  with  high  horizontal  resolution  and  low  vertical 
resolution,  a  Class  7  computer  (~  1  gigaflop  and  32  to  128  million 
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words)  is  required.  The  introduction  of  Class  7  computers  is 
anticipated  in  the  middle  and  late  1980s. 

Currently,  satellites  and  atmospheric  forecast  models  provide 
the  only  real  prospects  for  oceanographic  data  and  forcing  func- 
tions  with  global  coverage  and  resolution  adequate  for  use  in 
ocean  circulation  models.  For  this  purpose,  the  two  fields  of 
greatest  value,  which  are  observable  from  satellites,  are  alti¬ 
meter-derived  sea-surface  elevations  and  scatterometer-derived 
wind  stress.  Successful  ocean  prediction  appears  feasible  with¬ 
out  extensive  subsurface  data  acquisition,  by  using  the  simula¬ 
tion  skill  of  the  model  to  convert  the  potentially  well-observed 
fields  at  the  surface  into  subsurface  information.  This  approach 
looks  more  promising  than  attempting  to  assimilate  the  ex¬ 
tremely  sparse  subsurface  data,  a  process  likely  to  do  more  harm 
than  good.  Where  eddy-resolving  subsurface  data  with  useful 
coverage  are  available  (nowhere  at  present),  they  should,  of 
course,  be  used,  especially  if  they  are  available  on  a  regular  basis 
and  their  statistical  properties  are  well  known.  In  the  long  term, 
acoustic  tomography  shows  promise  in  this  area,  but  we  antic¬ 
ipate  that  subsurface  data  acquisition  adequate  for  prediction  of 
mesoscale  features  will,  at  best,  be  limited  in  coverage  for  at 
least  the  next  decade.  Some  subsurface  data  and  major  field  pro¬ 
grams  are  essential  for  the  local  and  regional  subsurface  vali¬ 
dation  of  the  forecast  models. 

The  surface-to-subsurface  data  conversion  process  would  be 
accomplished  by  using  both  the  altimeter  data  and  the  atmo¬ 
spheric  data  as  forcing  functions.  Although  this  is  an  untested 
hypothesis  at  this  time  and  the  assimilation  of  the  altimeter  data 
could  be  troublesome,  the  potential  for  surface-to-subsurface 
conversion  is  enhanced  by  noting  that  (1)  a  single  internal  vertical 
mode  explains  much  of  the  oceanic  variance;  (2)  in  a  layered 
circulation  model  with  a  single  internal  mode  there  is  a  one-to- 
one  correspondence  between  variations  in  the  sea-surface  ele¬ 
vation  and  variations  in  the  depth  of  the  pycnocline;  (3)  such  a 
model  can  be  initialized  by  altimeter  data  alone,  except  near  the 
equator,  where  wind  stress  is  necessary  for  initialization  in  ad¬ 
dition  to  forcing;  and  (4)  such  models  have  demonstrated  a  re- 
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markable  simulation  skill  in  certain  situations,  although  they 
miss  the  sometimes  important  effects  of  topography  and  baro- 
clinic  instability.  These  effects  require  at  least  one  additional 
vertical  mode;  we  recommend  no  less  than  three  in  an  opera¬ 
tional  forecast  model.  As  discussed  in  the  section  on  cost-ef¬ 
fective  design,  high  vertical  resolution  can  be  obtained  by  cou¬ 
pling  the  circulation  model  with  low  vertical  resolution  to  an  N 
x  1-D  mixed-layer  model  with  high  vertical  resolution. 

Altimeter  data  are  essential  for  global  prediction  and  repre¬ 
sentation  of  individual  current  meanders  and  eddies.  Approxi¬ 
mately  5-cm  accuracy  is  needed  for  variations  in  the  sea-surface 
elevation  related  to  ocean  surface  currents.  Track  spacing  of  not 
more  than  50  km  is  required  at  midlatitudes,  but  repeat  intervals 
up  to  a  month  could  probably  be  tolerated.  To  minimize  the  time 
required  to  span  an  eddy,  the  orbit  should  be  chosen  so  that 
adjacent  tracks  are  one  day  apart.  A  single  satellite  with  a  mul¬ 
tibeam  altimeter  and  a  20-day  repeat  orbit  could  meet  these  re¬ 
quirements  with  ease  (see  Bush  et  al.,  1984,  in  this  issue).  Min¬ 
imal  satisfaction  of  the  requirements  could  be  obtained  by  using 
a  single  satellite  with  a  nadir  beam  altimeter  and  a  40-day  repeat 
orbit  at  —1000-km  altitude,  if  a  numerical  prediction  model  is 
used  to  fill  in  the  temporal  gaps. 

In  the  absence  of  altimeter  data,  a  further  burden  would  be 
placed  on  the  simulation  skill  of  the  model.  Although  some 
oceanic  phenomena  are  quite  directly  related  to  atmospheric  or 
other  forcing  functions,  many  current  meanders  and  eddies  are 
not  because  they  can  result  from  flow  instabilities  and  because 
they  can  persist  for  many  months.  These  phenomena  could  only 
be  simulated  in  a  statistical  sense.  That  is,  the  circulation  model 
would  be  limited  to  forecasting  ensembles  of  eddies  and  the  gen¬ 
eral  character,  evolution,  and  movement  of  eddies  and  current 
meanders  in  different  regions. 

In  the  course  of  this  discussion  a  number  of  ideas  with  po¬ 
tential  application  to  ocean  prediction  have  been  encountered. 
Some  ideas  that  show  promise  are  (1)  the  use  of  satellite  altimeter 
and  scatterometer  data,  (2)  the  conversion  of  surface  into  sub¬ 
surface  information,  (3)  the  partial  decoupling  of  the  ocean  cir¬ 
culation  problem  from  the  turbulent  thermodynamic  mixed-layer 
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problem,  (4)  eddy-resolving  circulation  models  with  low  vertical 
resolution  and  basin-  or  global-scale  domains,  and  (S)  Class  7 
computers.  Some  that  are  not  promising  (at  least  within  the  next 
decade  or  more)  are  (1)  eddy-resolving  models  for  major  ocean 
basins  with  many  fixed  levels  in  the  vertical,  (2)  the  inclusion  of 
sophisticated  mixed-layer  physics  in  eddy-resolving  circulation 
models  covering  m^jor  ocean  basins,  (3)  the  assimilation  of 
sparse  in-situ  data  by  ocean  circulation  models,  (4)  coarse  grid, 
non-eddy-resolving  prediction  models  for  the  ocean  circulation, 
and  (5)  sub-basin-scale  models  with  extensive  open  boundaries, 
except  when  they  can  obtain  boundary  conditions  from  an  eddy¬ 
resolving  basin-scale  model.  Nowcasts  and  short-range  forecasts 
(a  few  days)  from  stand-alone  limited-area  models  appear  fea¬ 
sible  in  regions  with  adequate  subsurface  sampling  (none  at  pres¬ 
ent).  The  first  two  ideas  are  judged  not  to  be  promising  based 
only  on  anticipated  computing  power  within  the  next  decade  or 
more. 

Within  the  next  decade  the  appropriate  satellite  data.  Class  7 
computers,  and  eddy-resolving  basin-  or  global-scale  models 
should  become  available,  and  they  should  be  used  to  form  the 
heart  of  a  global  ocean  prediction  system.  In  the  meantime,  al¬ 
ternate  methods  described  in  this  paper  and  elsewhere  can  pro¬ 
vide  a  more  limited  capability  that  is  applicable  in  some  situa¬ 
tions.  In  particular,  forcing  functions  from  atmospheric  models 
should  permit  useful  nowcasting  and  predictive  skill  for  Class  1 
and  Class  3  phenomena  listed  in  Table  1  prior  to  the  availability 
of  satellite  altimeter  and  scatterometer  data. 
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